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The real-axis direct solution of the Eliashberg equations for the retarded electron-boson interaction in the 
half-filling case and in the presence of impurities is obtained for six different symmetries of the order parameter: 
s, s + id, s + d, d, anisotropics and extended-s. The spectral function is assumed to contain an isotropic part 
o% s F (fi) and an anisotropic one o? an F (Q) such that af 3 F (fi) = g ■ a\ n F (Q.) , where g is a constant, and the 
Coulomb pseudopotential ^t* is set to zero for simplicity. The density of states is calculated for each symmetry 
at T = 2, 4, 40 and 80 K. The resulting curves are compared to those obtained by analytical continuation of the 
imaginary-axis solution of the Eliashberg equations and to the experimental tunneling curves of optimally-doped 
Bi 2212 crystals. 



In this paper, we make use of the Migdal- 
Eliashberg theory jj] for the strong electron- 
boson coupling to discuss the effect of different 
possible symmetries of the order parameter [|| 
on the tunneling curves of copper-oxide super- 
conductors. Because of the layered structure of 
these materials, we can suppose the quasiparticle 
wavevectors k and k' to lie in the CuC>2 plane and 
call <f) and <f>' their azimuthal angles in this plane. 
Then we solve the Eliashberg equations (EE) us- 
ing a single-band approximation with a nearly- 
circular Fermi line. In the real-axis formalism 
the EE take the form of a set of coupled integral 
equations for the order parameter A(u),(f>) and 
the renormalization function Z(u>,<f)), containing 
the retarded interaction cv 2 (Sl, (f>, <f)')F(Q) and the 
Coulomb pseudopotential fj,* ((f), (f>') 

We hypothesize that the two last quanti- 
ties contain an isotropic and an anisotropic 
part ]3|,|j and we expand both of them in 
terms of basis functions. Actually, even 
though we are able to solve the EE for an 
arbitrary constant value of p* , we put it to 
zero for simplicity. The spectral function ex- 
panded at the lowest order is then expressed 
by: a 2 (Jl, (f), <f/)F(Sl) = a? s F(fi)^ s (<f>) ^ (</>') 
+ al n F(Cl)ip an {(f) ipan ((/>') where the basis 
functions tpi s ((f)) and tp an ((f)) are chosen as 
follows: ip is (<p)=l; ip an ((f>) =V2 cos (2(f>) for 
the d-wave, ifj an ((f)) =8-^/2/35 cos 4 (2(f)) for the 



anisotropic- s, and ip an (<fi)=— 2yj2/'d cos 2 (2(f>) 
for the extended-s 0. For simplicity again, we 
suppose that a 2 n F(n)=g-af s F(H,) where g is a 
constant ajm . Thus, the electron-boson coupling 
constants for the isotropic-wetve channel and the 
anisotropic-w&ve one, which are given by A<j S)an = 



(iM C d^i^) f + °° 



dn< arl F(tt)/ft, re- 



suit to be proportional: X an =9'\ 

We are interested in solutions of the real-axis 
EE of the form: A(w, (f))=A is (uj)+A an (uj)tp an ((f)) 
and Z(w,4>)=Z is (ui)+Z an (w)ip arl (<j>). The equa- 
tions are reported explicitly elsewhere Q. Here 
we suppose Z an (uj) to be identically zero [||J^]. 

The numerical solution of the real-axis EE is 
performed by using an iterative procedure. In 
view of a comparison to the experimental tunnel- 
ing curves obtained in Bi 2212 break junctions, 
we take a 2 s F(fi) = (A JS /A B i 2 2i2)a 2 ^(^)Bi22i2 
where a 2 F '(0)Bi22i2 has been experi- 
mentally determined in a previous paper ](| 
and rescaled to have T c =97 K. Once de- 
termined A(uj,(f>) and Z(tu,(f>), we calculate 
the quasiparticle density of states N(uj)= 
(1/2tt) J 27r d(f>Re (lu/^Jlu 2 - A(w,</>) 2 ) , whose 
convolution integral with the Fermi distribution 
is the quantity that must be compared to the 
experimental tunneling data. Incidentally, the 
explicit solution shows that the symmetry of 
A(u),(f>) is affected by the choice of the coupling 
constants Xi S and X an and, for some particular 
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Figure 1. Theoretical tunneling density of states for 
various symmetries at T=4 K and fi*=0. The circles 
represent our experimental data in Bi 2212 [[|. 

values of Xi S and A an , by the starting values of 
(uj) and A an (u). 
The figure reports the theoretical normalized 
conductance at T=4 K in the six symmetries an- 
alyzed and the Bi 2212 experimental data ob- 
tained in break-junction tunneling experiments 
||[|. The d-wave curve gives the best fit of the 
peak at the gap edge and of the conductance be- 
haviour inside the gap, suggesting that the sym- 
metry of Bi 2212 could be pure d-wave. Never- 
theless, none of the symmetries studied here is 
able to give, at the same time, T c =92-94 K and a 
peak at 35-41.5 meV, as experimentally observed 
by STM on Bi 2212 and recently reported in lit- 
erature (7). 

The same theoretical curves have also been ob- 
tained at T=2, 40 and 80 K. The anisotropic-s 
and extended-s curves are highly temperature- 
dependent, and some fine structures evidenced 



at T=2 K are already indistinguishable at 4 K. 
For T > T c /3 it is practically impossible to dis- 
tinguish between d-wave and mixed symmetries, 
while the s-wave curve remains clearly distinct. 

The curves at T=2 K can be compared to 
those obtained by analytically continuing the 
imaginary-axis solutions A(iw„) and Z(iw n ). In 
general, at this low temperature, the analytical 
continuation gives a reasonable agreement with 
the real-axis solutions. The agreement is satisfac- 
tory for the s + id, d, anisotropic-s cases, and be- 
comes a little worse in the extended-s case. The 
imaginary- axis s-wave curve is instead markedly 
shifted (of about 3 meV) toward higher energies, 
and then is unable to approximate the real-axis 
solution. In the (s + d)-wave symmetry, even the 
shape of the curve is heavily different in the two 
cases B. 

Finally, we have studied the effect of 
non-magnetic impurities on the tunneling den- 
sity of states by solving the appropriate real-axis 
EE We have found that this effect is greater 
on the d-wave component of the order parame- 
ter. A small amount of impurities in the uni- 
tary limit gives rise to a zero bias in the d-curve, 
and then improves its agreement with our ex- 
perimental points |J. In the non-unitary case, 
the d-curve peak is further lowered, broadened 
and shifted leftward. In the same conditions, all 
the mixed-symmetry curves are modified in non- 
trivial ways, e.g. the s + id and the anisotropic-s 
curves become practically indistinguishable from 
the s-wave one, but are shifted toward lower en- 
ergies. Finally, the s-wave curve is nearly unaf- 
fected by this kind of impurities, apart from a 
small shift of the gap toward higher energies. 
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